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ABSTRACT 

We have compared the B — V colour excess, E(B-V), obtained for a sample of five opti- 
cally visible massive YSOs both from diffuse interstellar bands (DIBs) in their spectra 
and from their optical continuum slopes. Our targets are HD 200775, BD+40°4124, 
MWC 1080, MWC 297 and MWC 349A. First, E(B-V) towards each of the targets is 
derived by dereddening the observed continua to match those of B-type standard stars. 
A survey of DIBs in the spectra of the massive YSOs, and a control field star, then 
reveals that the DIBs are significantly weaker in the former than would be expected 
based on the total E(B-V) values. This result is strengthened by the finding that the 
DIBs in the control field star, HD 154445, have on average the strength expected from 
its continuum E(B-V). 

A rough estimate of the foreground reddening of intervening diffuse interstellar 
medium shows it to be smaller than the DIB- E(B-V), implying that at least part 
of the DIB carriers are formed within the parental molecular clouds in which the 
YSOs are embedded. The formation efficiency of the DIBs varies strongly however 
from cloud to cloud. The DIB- E(B-V) compares favourably with the total E(B-V) 
towards BD+40°4124, but is almost negligible in the line of sight towards MWC 297. 
Despite this general, but not unexpected, deficit we provide evidence that the DIB at 
5849 A is a good tracer of total extinction in these lines of sight. 



Key words: stars: circumstellar matter - stars: emission line, Be 
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1 INTRODUCTION 

Still unidentified, diffuse interstellar bands (DIBs) are be- 
lieved to be the result of the absorption of starlight by 
complex molecules or small dust particles in the interstellar 
medium. The weak absorption features have been well doc- 
umented by early spectroscopists although it was not until 
the 1930s that the interstellar nature of the lines was recog- 
nised by Merrill (1934) and coworkers. Continuing investiga- 
tions showed that the strength of the interstellar absorption 
features increased with the star's distance and degree of in- 
terstellar reddening (Merrill, 1936) . Hence the measurement 
of the strength of the DIBs is related to path length through 
the diffuse interstellar medium (ISM). More recently it has 
become clear that the strength of the DIBs correlate with 
the column density of Hi, but not with the H2 column (Her- 
big 1993, 1995). 

Although the line strengths of individual DIBs correlate 
relatively well with the foreground reddening in cases where 
the diffuse ISM is the main contributor to the extinction, 
there are particular classes of sightline where this correlation 
breaks down. For example, for lines of sight through circum- 



stellar material the DIB-derived colour excesses ( E(B-VJ) 
turn out to be smaller than the total E(B-V) determined 
by other means (e.g. Snow & Wallerstein, 1972; Le Bertre 
& Lequeux, 1993 on mass-losing objects; Porceddu, Ben- 
venuti & Krelowski 1992, on Be stars). A similar effect has 
been found for sightlines into star-forming regions, with the 
further complication that not all DIBs suffer the same re- 
duction (e.g. Meyer & Ulrich 1984, who studied T Tau ob- 
jects; Adamson et al. 1991, Taurus dark clouds; Jenniskens, 
Ehrenfreund & Foing 1994, Orion). These findings suggest 
that the physical conditions in these media are such that the 
carriers responsible for the DIB absorption do not exist in 
the same relative quantities as in the diffuse ISM, which is 
probably related to the formation, excitation and destruc- 
tion of different DIB carriers. The concept of different DIB 
'families' (Krelowski & Walker, 1987) can be traced back to 
this concept. Some band strengths are well correlated with 
each other, but not with the reddening, while other bands 
correlate strongly with reddening. 

In the present study we describe the behaviour of the 
DIBs in the lines of sight towards massive young stellar ob- 
jects. Our targets are the optically visible young stellar ob- 
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ject (YSO) MWC 349A and four Herbig Be stars, objects 
which are optically bright enough to obtain high signal-to- 
noise spectra at intermediate dispersion in the optical. The 
Herbig Ae/Be stars are widely believed to be associated with 
circumstellar envelopes of dust. The material in the line of 
sight towards these young massive objects is thus a combi- 
nation of 'normal' intervening diffuse interstellar material, 
possibly cold material in the clouds in which the objects 
have formed, and certainly warm circumstellar material. It 
is hard to estimate the relative contributions to the total 
extinction of these media, but observations of the DIBs can 
in the long run help to elucidate this problem. 

Here, we investigate whether the DIBs in the spectra of 
our targets would imply the same extinction as that derived 
by other means. Any difference between the derived E(B- 
V) values would set an upper limit to the contribution of 
diffuse material to the line-of-sight extinction. In the follow- 
ing we will derive E(B-V) by comparing the observed slopes 
in the target spectra with standard star spectra - a 'contin- 
uum method', and derive E(B-V) from the DIB strengths 
alone using published relations between the EW of DIBs and 
E(B-V). 

The lay-out of this paper is as follows: In SecJd we de- 
scribe the observations and data reduction. In Sec. ^ the 
determination of the continuum reddening towards the tar- 
gets is presented and compared with values in the literature. 
In Sec. ^| the measurement and analysis of the DIBs in the 
spectra are presented. We then compare the DIB derived 
E(B-V) with statistical estimates of the colour excess at- 
tributable to the foreground ISM described in Sec. [B|. The 
paper ends with a summary. 



2 OBSERVATIONS AND DATA REDUCTION 

In the course of our study of the optical spectra of mas- 
sive YSOs (e.g. Drew et al. 1997) we obtained observations 
on the nights of the 20th to the 22nd June 1994 with the 
William Herschel Telescope (WHT) at the La Palma Obser- 
vatory. The weather throughout the run was good, yielding 
a seeing of 0.7 to 1 arcsec. The twin-beam intermediate dis- 
persion spectrograph (ISIS) was used with the R1200B and 
R1200R gratings on the blue and red arms of the instrument. 
The blue detector was a 1124 x 1124 pixel TEK1 chip, the 
red detector was a 1180 x 1280 EEV6 chip. The wavelength 
ranges covered were 3860-5190A in the blue arm with a 
small gap of 120A between 4270-4390A . The red arm cov- 
ered the wavelength range 5840-8390 A. The dichroic crossed 
over at 5400 A. The instrumental set-up does not allow for 
the observation of the well-known DIBs at 5780 and 5797 
A. However, the large wavelength coverage and high sensi- 
tivity of the observations make up for this loss. For wave- 
length calibration a copper-neon-argon lamp was observed 
at frequent intervals. The projected entrance slit width was 
1 arcsec, which resulted in spectral resolution elements of 
0.8A (A/ A A ~ 8200 at Ha), as determined from arc line 
profile fits. To minimise cosmic ray events in the spectra, 
exposure times were limited to 1800s. The observations in- 
cluded a spread of spectral type standards which are listed in 
Table |l|. The typical integration times for these bright stars 
were of order minutes per setting, resulting in SNR ratios 
in the continuum of more than 100. Table \A lists the YSOs 



Table 1. The spectral type standards and atmospheric standards. 



Object 



SpTp 



HD 154445 
HD 186618 
HD 214432 
HR 6092 
HR 8335 



B1.5V 

B0.7IV 

B2.5V((n)) 

B5IV 

B3III 



The spectral types are taken from Walborn (1971) 

observed. At each wavelength setting the spectral coverage 
obtained was generally a few Angstroms in excess of 400A. 
The SNR achieved for each of the objects is given in the 
final column in the tables. 

Data reduction was performed using iraf. Additional 
manipulation made use of both the iraf and STARLINK- 
dipso software packages. The data were reduced from two- 
dimensional frames into a one-dimensional form, through the 
steps of de-biasing, flat- fielding, sky subtraction and wave- 
length calibration. A first order cubic spline was used for 
the wavelength calibration. The fit attained an rms error 
in the residuals of less than 0.03A. The extracted spectra 
were airmass corrected in order to remove the effects of at- 
mospheric extinction in the spectra. We applied the appro- 
priately scaled sum of the wavelength-dependent mean ex- 
tinction curve for La Palma (King 1985) and the "grey" 
aerosol component whose magnitude is measured nightly 
and archived by the Carlsberg Meridian Circle. Where nec- 
essary, telluric features were removed by dividing the tar- 
get star spectra by the corresponding spectra of a suitable 
comparison star (either HR 6092 or HR 8335). The spec- 
tra were scaled such that the depths of the telluric absorp- 
tion features in the standards match those in the target 
spectra. The main regions affected by telluric features are 
5860-6000A, 6270-6330A 6450-6580A and 6860A-6950A. 
All wavelengths longwards of 6950A are susceptible to spo- 
radic atmospheric features. 



3 CONTINUUM DETERMINATION OF 
COLOUR EXCESS. 

Before assessing the relation between the DIB strengths and 
reddening towards the target stars, we need to have an esti- 
mate of the total E(B-V) for each of them. In this section we 
will derive E(B-V) by comparing the observed continuum 
slopes of the spectra with those of spectral standard stars. 
Their spectra were observed during the same run with the 
same instrumental set-up. All observations were obtained 
with the slit at a parallactic angle to minimise the influence 
of the dispersion of the Earth atmosphere. As stated before, 
the spectra are corrected for airmass. The procedure is thus 
independent of the observed flux from the objects, and only 
takes into account changes in the spectral slopes. It is the 
same as that used by Cohen et al. (1985) and Drew et al. 
(1997) to derive the reddening to MWC 349A and MWC 
297 respectively, both stars are also included in this sample. 

The extracted spectra of the target stars are dered- 
dened to find the best match to the continuum slopes of 
the standard stars, with E(B-V) as the only free parame- 
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Table 2. Observations of the optically bright YSOs. The effective 
wavelength coverage was approximately 400A per setting centred 
on the wavelengths given in the second column. The co-added 
SNR is given if more than one observation is available. The B 
and V magnitudes are taken from Hillenbrand et al. (1992), the 
V magnitude for MWC 349A is taken from Cohen et al. (1985). 



Table 3. Determination of the total E(B-V). The central wave- 
length of the spectral range is given in column 2. Column 3 records 
the values of the E(B-V) as determined by the H83 law derived 
using each of the standard stars. Below the lines of measured 
E(B—V)b the weighted mean of all the E(B-V)s is presented. 
Previous E(B-V) or Ay obtained from the literature are given 
in the final column. 



Object 
Spectral type 



V Wavelength 
(centre, A) 



Exposure 
Time (s) 



Continuum 
SNR= 



HD 200775 
B2.5e 4 



7.75, 7.42 



BD+40°4124 
B3C 1 



11.28, 10.54 



MWC 1080 
BO 2 



13.08, 11.68 



MWC 297 
B1.5Ve 3 



14.54, 12.17 



MWC 349A 



13.5 



4600 
5000 
6055 
6475 
6895 
7315 
7735 

4600 
5000 
6055 
6475 
6895 
7315 

4600 
5000 
6055 
6475 
6895 
7315 

4600 
5000 
6055 
6475 
6895 
7315 
7735 

4600 
6055 
6475 
6895 
7315 



50 
50 

100 
30 

100 
70 
10 

1200 
600 
600 
630 
200 
200 

2648 
1000 
1800 
200 
300 
300 

9000 
1800 
5400 
2530 
2100 
1200 
500 

5400 
800 
270 
1406 
2500 



Object 



E(B-V) 
B0.7,B1.5,B2.5 



Literature 



Method 



239- 
250 
340 
140 
390 
300 
220" 

210 
130 



HD 200775 



4600A 
5000A 



0.60, 
0.75, 



0.45, 0.55 
0.45, 0.70 



mean: 0.6 



0.65 1 
0.56 2 
0.64 3 
0.57 4 



colours 
colours 
SED model, 
colours 



BD+40°4124 



4600A 
5000A 



1.00, 0.85, 0.95 
1.15, 0.90, 1.00 
mean: 1.0 



0.96 1 
0.94 2 



colours 
colours 



190 MWC 1080 
150 
120 
135- 



4600A 
5000A 



1.60, 1.45, 1.45 1.71 1 
1.80, 1.50, 1.75 A^=5.42 
mean: 1.6 



colours 
as here 



110 
100 
200" 
40 
110 
135 



MWC 297 4600A 2.95, 2.85, 2.95 2.68 1 
5000A 3.05, 2.70, 2.85 A v =7.9 ' 
mean: 2.9 



colours 
NIR slope 



MWC 349A 4600A 



2.95, 2.95, 3.00 A v =9.9 7 as here 
mean: 3.0 A v =8.8 6 NIR slope 



150 



References: 1 - Hillenbrand et al. 1992, 2 - Finkenzeller & 
g( Mundt 1984, 3 - Voshchinnikov, Molster & The 1996 4 - Pfau et 
al. 1987, 5 - Cohen & Kuhi 1979, 6 - McGregor, Persson & 
Cohen 1984. 7 - Cohen et al. 1985 



References to the spectroscopic determinations of the spectral 
types: 

1 : Finkenzeller (1985), 2 : Cohen & Kuhi (1979), 3 : Drew et al. 
1997. 

ter. This method implicitly assumes that the intrinsic con- 
tinuum spectral energy distributions of the Herbig Be stars 
in the wavelength ranges of interest resemble those of the 
spectral standards closely. This is indeed the case for our 
targets. Spectroscopic determinations of their spectral types 
(Table ^) cover the same spectral type range as our compar- 
ison stars (see Table [l]). We are not aware of a spectroscopic 
spectral type for the extreme emission line star MWC 349A 
- but considering its emission line character and continuum 
energy distribution, it is widely held to be an early B type 
star (see e.g. Cohen et al. 1985). The knowledge of the pre- 
cise spectral type is not critical however, since tests with 
unreddened Kurucz (1991) models indicate that the differ- 



340 
300 
320 
305 
300 

^.ence in slopes between the extreme spectral types of 09 and 
120^ corresponds to an error in the derived E(B-V) of be- 
-Ql^ween 0.05 and 0.1, which is much less than our adopted 
2Qierror in the final E(B-V) values. 

355 Additionally, account needs to be taken of the fact that 
the standard stars are lightly reddened themselves. We de- 
riv ed their reddening using the observed colours (obtained 
from the SIMBAD database) and intrinsic colours listed in the 
tables of Schmidt-Kaler (1982) for their respective spectral 
types. HD 154445 has an E(B-V)=0A1, while HD 186618 
and HD 214432 display a small amount of reddening of 0.07 
and 0.11 respectively. Dereddening HD 154445 with respect 
to HD 186618 and HD 214432, using the same technique as 
for the young stellar objects, yielded an E(B-V) difference 
consistent with the photometric estimates. 

Since we are only interested in the continuum slope, 
strong emission lines (in the Herbig Be stars) and absorption 
lines (in the standard stars) were snipped out of the spectra 
to prevent these from affecting the fitting procedure. 

Once the spectra had been prepared, the fit proce- 
dure was relatively straightforward: each spectrum was de- 
reddened for various values of E(B-V) until an rms mini- 
mum between the de-reddened spectrum and the standard 
star was reached. The procedure encountered too few points 
upon which to operate satisfactorily in the case of extreme 
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emission line stars MWC 1080 and MWC349A. In these 
cases, the colour excesses were better determined by eye. 

For the determination of the reddening we have used 
reddening laws from two sources, the mean Rv = 3.1 galactic 
extinction law compiled by Howarth (1983, H83), and also 
the Rv dependent analytic prescriptions due to Cardelli, 
Clayton & Mathis (1989, CCM). The red settings proved 
not sensitive enough to return useful extinction estimates, 
and only the results for the 4600 and 5000 A settings are 
listed in Table |. 



3.1 Results 

The YSO colour excesses, derived using the H83 law, are set 
out in Table ^. These include the correction for reddening of 
the standard stars. Below the lines of measured E(B-V)s, 
the weighted mean of all the E(B-V)s for a given star is 
presented. As we associate 3<r errors of 0.15 and 0.2 with the 
4600 and 5000 A settings respectively, the weights applied 
to them in obtaining the mean were chosen to be in the 
ratio 4:3. It is worthy of note that the derived E(B-V) of 
MWC 349A and MWC 1080 are in good agreement with the 
literature values despite the extreme emission line character 
of their spectra. 

In principle the CCM law follows - as the H83 law - 
the van de Hulst no. 15 law (van de Hulst, 1949) closely 
for Ry= 3.1. The main advantage of the CCM law over 
the H83 law is that the dependence on the ratio of total 
to selective reddening may be investigated. It turned out 
that the Rv sensitivity in our spectra is very small; fits 
using different values for Rv yielded similar fit qualities and 
results as those with Rv = 3.1. On the whole, the CCM law 
results were within la of the results outlined in the table. 
Exceptions to this were obtained for the highly reddened 
objects MWC 349A and MWC 297 where the CCM law 
yielded lower E(B-V) values by as much as 0.45 in the 
4600 A setting. This is because the CCM law is somewhat 
steeper in the 4600 A wavelength range than the H83 law, 
a difference that is amplified for highly reddened objects. 



4 DIB DETERMINATION OF COLOUR 
EXCESS 

We selected all DIBs listed by Herbig (1995) which appear in 
the wavelength range covered by our data. The DIBs in Her- 
big's original list with a normalised equivalent width < 8 mA 
are below our detection capabilities and are not considered. 
To this list we added stronger lines listed by Jenniskens & 
Desert (1994, hereafter JD) but not by Herbig (1995). The 
final list comprises 80 DIBs of which 13 are from the JD 
catalogue. 

With the search list in hand, the spectra were ex- 
amined independently for DIBs by two of us to minimize 
the inevitable subjective element. In general, the two in- 
dependent assessments yielded identifications and equiva- 
lent widths that were in agreement at the 10-15% level. In 
those cases where extreme differences appeared, such as non- 
identifications or very different equivalent widths, the spec- 
tra were inspected again. Often, the combination of faint fea- 
tures and the presence of emission lines in the spectral neigh- 



1.8 



1.6 



1.4 



MWC 349A 




-| 2 - HD 200775 



- HD 154445 



J_ 



J_ 



5840 5845 5850 5855 

Wavelength (A) 



1 1 1 1 I 1 1 1 'y^ 1 

Atm. Standard 



1.5 




HD 154445 



J_ 



J_ 



7215 7220 7225 7230 

Wavelength (A) 



Figure 1. Continuum normalized spectra around the A5849 and 
A7224 DIBs. The stars are plotted in order of decreasing E(B- 
V) from top to bottom. The reddened standard star HD 154445 
is shown as comparison in both figures. For the A7224 DIB, the 
atmospheric standard HR 8335 is shown to illustrate results of 
the correction of telluric absorption lines. 
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bourhood caused the differences. The measured FWHM are 
in accordance with their listed values. 

Because of the emission line nature of the objects, many 
lines were hampered by blends with emission lines. The spec- 
tral resolution of the observations is in most cases enough to 
sample the emission lines sufficiently to enable us to distin- 
guish them from the DIBs. After carefully selecting the least 
affected lines, a sample of 18 DIBs remained. These lines are 
presented in Table ^. The reddened spectral type standard 
HD 154445 displays a number of DIBs and is also included 
in the investigation as a fortuitous control. Examples of two 
DIBs are given in Fig. 1. 

In the following, we will discuss the reddening as traced 
by the DIBs, and exploit the concept of a DIB- E(B-V). 
These are calculated using the conversions EW- E(B-V) 
from JD. Although often the EW- E(B-V) calibration of 
Herbig (1995) is used, JD provide the mean normalization 
over four different lines of sight sampling just the diffuse 
ISM, instead of only one object (HD 183143), which reduces 
the effect of possible peculiarities in using one particular line 
of sight. While EW/ E(B-V) for individual DIBs can vary 
factors of several from sightline to sightline in the data of JD, 
these variations smooth out to a scatter of 10-20% around 
the mean when a larger sample of DIBs, such as ours, is 
considered. This suggests that calculation of E(B-V) from 
a sample of DIBs gives a good indication of the minimum 
contribution to the total reddening along the sightline by 
the diffuse medium (cf. Herbig 1995). 

We note in passing that Herbig's EW values for the 
DIBs in HD 183143 are larger than those listed by JD. In 
fact, JD also note that the EW values of HD 183143 pre- 
sented by Herbig (1975) and Herbig & Leka (1991) 'may' be 
systematically larger by 20%. Le Bertre & Lequeux (1993) 
also mention that their measurements of the DIBs of HD 
183143 'possibly give smaller equivalent widths' than Her- 
big. We find the same in our own high resolution UES spec- 
trum of HD 183143 (courtesy of Ton Schoenmaker and Eric 
Bakker) . 

4.1 Results 

The results are summarized in Fig. ^, where the inferred 
E(B-V) of the DIBs is plotted against the continuum E(B- 
V) listed in Table || Note the panel showing the relation for 
the A5849 line. It is clear that this line results in deduced 
E(B-V) values for every star (except for HD 200775 where 
the line is not detected) which are in agreement with the 
continuum-determined E(B-V). The spectra around this 
line are shown in Fig. |l| Our data strengthen the tentative 
suggestion of Chlewicki et al. (1986) that the A5849 line 
is frequently a good indicator for the total extinction to- 
wards an object. A crude match between line- E(B-V) and 
continuum- E(B-V) also seems to exist for the AA4963, 4726 
bands. However, in most other lines, the trend is that the 
line- E(B-V) value is lower than that determined from the 
continuum. 

An illustrative plot is provided in Fig. || It shows the 
inferred E(B-V) of the measured lines against DIB wave- 
length. The horizontally drawn lines indicate the continuum 
E(B-V) (Table |). In especially BD+40°4124 and the stan- 
dard star HD 154445, the DIBs return an E(B-V) close 
to the value of the total E(B-V), but in the other cases it 



Table 5. E(B-V) as derived from the DIBs. n is the number of 
lines used in the calculation 





Cont. 




DIB 


Star 


E(B-V) 


n 


< E(B-V)> 


HD 154445 


0.4 


15 


0.43 ( 0.07) 


HD 200775 


0.6 


9 


0.14 ( 0.07) 


BD+40°4124 


1.0 


18 


0.64 ( 0.08) 


MWC 1080 


1.6 


17 


1.01 ( 0.10) 


MWC 297 


2.9 


18 


0.47 ( 0.08) 


MWC 349A 


3.0 


16 


1.56 ( 0.11) 



can be grasped straightaway that the DIB E(B-V) is sig- 
nificantly smaller than that determined for the continuum. 
The three DIBs in the spectrum of MWC 297 that are suffi- 
ciently strong to trace the total E(B-V) are the same three 
mentioned above as correlating well with the total E(B-V). 
All the other DIBs return values very much smaller than the 
continuum E(B-V). 

The mean DIB- E(B-V), calculated from all detected 
lines is presented in Table H. In order to avoid a bias to- 
wards the larger EWs, the mean has been calculated using 
a,l/<T weighting, instead of weighting the measurements ac- 
cording to the inverse square of their errors (the officially 
commended weighting scheme). We did not take into ac- 
count any non-detections, but it can be seen in Fig. ^ that 
the resulting upper limits on inferred E(B-V) values are 
mostly in agreement with the means obtained. 

In the only instance where we are certain that the dif- 
fuse ISM is the principal extinction medium, HD 154445, 
the DIB E(B-V) is in pleasing agreement with the contin- 
uum E(B-V). For the massive YSOs, it appears that the 
DIB- E(B-V) is always smaller than the continuum E(B- 
V), but that the extent of the shortfall differs substantially 
from object to object: specifically, for BD+40°4124 the DIB- 
E(B-V) is about 65% of the continuum E(B-V), while the 
other extremes are MWC 297 and HD 200775, where the 
DIBs seem to trace only 20% or less of the total E(B-V). 



5 THE CONTRIBUTION OF FOREGROUND 
REDDENING TO THE TOTAL EXCESS 

Now that the total E(B-V) and the E(B-V) traced by the 
DIBs are derived, it is of interest to have an independent 
estimate of the extinction due only to the foreground diffuse 
interstellar medium along our sightlines. This is done by 
investigating the reddening of nearby field stars. 

To facilitate this, we have searched the polarization cat- 
alogue compiled by Matthewson et al. (1978), which con- 
tains, in combination with polarization data, values for the 
extinction and photometric distances of more than 7500 ob- 
jects. For each of our stars, we selected all objects within 
a radius of 400 arcmin. The results of this exercise are set 
out in Fig. ^, where the E(B-V) for each star is plotted 
against its distance modulus m-M. A consistency check was 
made by a comparison with similar relations presented by 
FitzGerald (1968) and Neckel & Klare (1980), and by calcu- 
lating the E(B-V) and distance modulus using the spectral 
types, intrinsic B-V, and Mv (using Ry= 3.1) listed by 
Schmidt-Kaler (1982) for randomly selected objects. This 
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Table 4. The DIBs in the spectra of the 5 massive YSOs and the control star HD 154445. The first two columns identify the DIBs 
by their wavelength and their respective normalization to unit E(B-V) as published by JD. Errors on the equivalent width are given 
between brackets. The value for the A6203 DIB is the sum of the 6203 and 6204 A components resolved by JD, but not by us. NO: Not 
observed, B: blend with emission line, an idication of the EW is given NP: Not present — in general an upper limit of 10mA is found. 



A W/ E(B-V) HD 200775 HD 154445 BD+40°4124 MWC 1080 MWC 297 MWC 349A 

A mA 



4428 


2231 


330 ( 70) 


890 ( 80) 


1070 ( 270) 


1740 ( 450) 


1790 ( 360) 


3420 ( 480) 


4726 


123 


10 ( 5) 
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Figure 2. Deduced E(B-V) from the measured EWs of the lines for the different objects plotted per line. The open circles represent 
the bands affected by emission lines. The solid lines indicate the E(B-V) (DIB) equals the E(B-V) (continuum) relation. 



showed that the listed E(B-V) and rn-M in the Matthew- at their estimated distances |j and horizontal lines represent- 
son et al. (1978) catalogue are in general correct to within 
0.2 and 0.5 magnitudes respectively. 

* The adopted distances to the objects are as follows: MWC 
297 at 250 pc based on its spectral type (Drew et al. 1997); 
In the figure, the objects are indicated with vertical lines MWC 349A - at least 1000 pc, based on spectral type of com- 
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Figure 3. Inferred E(B-V) from the DIBs on a line by line 
basis. The horizontal lines represent the continuum E(B-V). 



ing the DIB- E(B-V). The lines of sight to BD+40°4124, 
MWC 349A and MWC 1080 initially show a trend between 
the distance modulus and E(B-V), which then breaks down 
at large distances, where the objects are scattered in E(B- 
V). Since all these objects are located in, or near, molecular 
clouds the scatter is real, and reflects the many stars at 
different depths and locations within patchy clouds. Where 
feasible, a conservative upper bound to the foreground red- 
dening in the lines of sight is indicated by the solid lines. 

It is not ruled out that the sightlines to MWC 297 and 
HD 200775 show a similar behaviour to that seen in e.g. 
MWC 349A, but the sparse number of data points makes it 
hard to recognize this. For example, MWC 297 is known to 
be located within the Aquila Rift at a distance of about 250 
pc (Dame & Thaddeus, 1985, Drew et al. 1997). However, 
one may see that it is only beyond ~500 parsec (m-M ~ 
8.5) that the data points scatter in a similar way as for 
MWC 349 A. These stars must be behind the Aquila Rift. 
It would appear that the catalogue of Matthewson et al. 
includes too few stars that are located in the Rift at m- 
M ~ 6-7. Similarly, it seems in the figure for HD 200775 
(positioned within NGC 7023) that E(B-V) only starts to 
become significant beyond m-M ~ 6. 

It is clear that in all cases the foreground extinction 
is much less than the total extinction towards the objects 
(Fig. [| Table ||). For the E(B-V) traced by the DIBs some 
interesting differences appear. For both MWC 297 and HD 
200775 the foreground extinction may be comparable with 
the E(B-V) as calculated from the DIBS, but not larger. In 



panion MWC 349B (Cohen et al. 1985) ; HD 200775, 600 pc by 
Rogers, Heyer & Dewdney (1995) based on spectral type B2.5IV. 
BD+40°4124, 1000 pc, Hillenbrand et al. 1995; MWC 1080, 2500 
pc, Canto et al. (1984), but is likely to be smaller - see Hillen- 
brand et al. (1982), who cite 1000 pc. HD 154445, 275 pc, based 
on spectral type. 
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Figure 4. The relationship between E(B-V) and distance mod- 
ulus of objects around the position of the target stars. The circles 
indicate objects within 200 arcmin, the crosses represent objects 
within 400 arcmin. The locations of the target stars are indicated 
by the intersection of their DIB- E(B-V) and distance. Note that 
the total E(B-V) to the objects is in all cases larger than the 
DIB- E(B-V) (Table ^).The solid lines represent a conservative 
estimate of the contribution of the foreground reddening as a 
function of distance. 



the other cases, the DIB- E(B-V) is intermediate between 
the total E(B-V) and the foreground extinction. This sug- 
gests that the molecular clouds, or even the circumstellar 
material, harbour conditions in which the DIB carriers can 
be present albeit at reduced abundance. 

We remind however that there is a small collection of 
DIBs in the spectrum of MWC 297 that is particularly 
strong; namely A5849, A4963 and A4726 appear to trace the 
total continuum E(B-V). 



6 CONCLUDING REMARKS 

DIB features are in general well-correlated with the inter- 
stellar extinction for field objects, but deviate in instances 
where circumstellar or parental cloud material is present. 
The results in Table || and in the previous section certainly 
indicate that the latter is true for our sample of massive 
YSOs as well. The diffuse ISM contributes only a fraction of 
the total reddening towards the YSOs - the rest of the ex- 
tinction being due to either cloud or circumstellar material. 

In conclusion, our study of Diffuse Interstellar Bands 
towards a sample of massive YSOs can be summarized as 
follows; 

(i) From our crude estimates of the foreground extinction 
we find that a large fraction of the extinction to these stars 
is not provided by the diffuse interstellar medium, but by 
circumstellar material or dust within the parental cloud. 

(ii) By and large, the DIB- E(B-V) is smaller than the 
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total E(B-V) towards the optically visible massive YSOs - 
but is mostly larger than the expected foreground redden- 
ing. This means that at least a fraction of the DIBs are 
formed within the molecular clouds, but less efficiently than 
in the normal diffuse ISM. This implies that conditions in 
this material are such that DIB excitation is weaker or that 
the DIB carriers are less efficiently formed or more prone 
to destruction. The DIBs trace most of the E(B-V) suf- 
fered by BD+40°4I24 but very little in the line of sight to 
MWC 297. Such variations in 'efficiency' clearly imply that 
the conditions vary strongly from one cloud of massive star 
formation to the other. 

(iii) If the DIBs are due to different types of carrier, the 
strength of the lines will be proportional to the respective 
column densities of the carriers. The effective normalization 
of line-strengths to E(B-V) for 'normal' diffuse interstellar 
material then gives us a handle on their relative abundances. 
From the fact that the A5849 line seems to trace the total 
E(B-V), one may infer that the carrier of this line is more 
readily formed or less likely to be destroyed, in star-forming 
clouds than other lines. The same may also be true for the 
A4726 and A4963 DIBs as strikingly illustrated by the case 
of MWC 297. 
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